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Abstract: 

Supersymmetry predicts that gauge couplings are equal to the corresponding gaugino- 
sfermion-fermion Yukawa couplings. This prediction can be tested for the QCD sector 
of the MSSM by studying the processes — > squark+antisquark+gluon and e + e~ — ► 
squark+antiquark+gluino at a future linear collider. We present results for these processes at 
next-to-leading order in a s in the framework of the MSSM. We find sizable SUSY-QCD correc- 
tions. The renormalization scale dependence is significantly reduced at next-to-leading order. 
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1 Introduction 



Softly broken supersymmetry predicts that each known particle has a superpartner with equal 
gauge quantum numbers and spin different by 1/2. Further, the gauge couplings have to be 
equal to the corresponding gaugino-sfermion-fermion Yukawa couplings. These coupling re- 
lations are vital for the cancellation of quadratic divergencies. Therefore, in order to establish 
supersymmetry experimentally, one not only has to find new particles and measure their quan- 
tum numbers, but also to verify the SUSY coupling relations. For the QCD sector of the MSSM, 
the coupling relation is depicted in Fig. 1. Analogous relations hold for the electroweak sector 
©• 

The SU(3) coupling relation can be verified at colliders by measuring the cross sections 
for the processes 
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and compare them to precise theoretical predictions, i.e. to a calculation at next-to-leading order 
(NLO) in oCj. The standard NLO QCD corrections to the 3-jet production process (§a) are well 
known O]. Here we will also present the virtual SUSY-QCD corrections to order oc^ for this 
process. Cross sections for processes (Qb) and (§c) have been computed to leading order in [Q]. 
We have recently performed a full NLO calculation for these processes. In this talk, we will 
discuss our main results. A detailed account of our work will be given elsewhere [Q]. 



2 Squark and gluino production at leading order 

We first discuss the leading order cross sections for processes (Qb) and (Qc). We allow only 
light quarks q = u,d,s,c,b (and antiquarks) in the final state of process (^c) and neglect their 
masses. We further exclude scalar top quarks as final state particles. The mixing between the 
chiral components qi and qn is neglected and all five squark flavours are assumed to be mass 
degenerate. The cross section for (Qb) diverges as the gluon energy goes to zero. We define an 
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Gauge coupling g s = Yukawa coupling g s 

Figure 1 : Coupling relation for the QCD sector of the MSSM. 
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infrared finite cross section as follows: 

o(E cut )= H o(e + e~ -> q h Z[ h X\ E x > E cxA ), (2) 

q=u,d,s,c,b h=L,R 

where X = g at leading order. The minimal energy E cut can be chosen by the experimentalist. 
We will take E cm = 50 GeV for most of the numerical results below. For process ([2|c), the total 
cross section is defined by 

<*tot= £ H jo(e + e~ ^ q h qg)+o(e + e~ ^ (3) 

q=u,d,s,c,b h=L,R 

Figs. 2a and 2b show the leading order cross sections o(£ , cut ) and a to t defined above. We also 
plot for comparison the LO cross section for e + e~ — > qq. We distinguish between the cases 
m q < nig (Fig. 2a) and m q > nig (Fig. 2b). One sees that even for rather large squark and gluino 
masses, the cross section o(£ , cut ) defined in @ reaches values of several tens of femtobarn, 
while the cross section a to t defined in (0) is about one order of magnitude smaller. 




3 SUSY-QCD corrections to e + e —> qqg and e + e ^qqg/qqg 

The calculation of the NLO SUSY-QCD corrections to processes (@b) and (0c) requires the eval- 
uation of both virtual and real corrections. Possible real corrections at order to process (Qb) 
are contributions from the processes e + e~ — > qqgg and e + e~ — > qqqq, while the real corrections 
to process (|2|c) consist of real gluon emission, e + e~ — > qqgg . Sample Feynman diagrams for 
both virtual and real corrections are shown in Figs. 3a-3d. 
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Figure 3: Sample Feynman diagrams for virtual and real corrections to e + e — > qqg (Figs, (a) 
and (b)), and to e + e~ — > qqg (Figs, (c) and (d)). The trivial leptonic part is not shown. 



In the calculation of the virtual corrections one encounters both ultraviolet and soft and collinear 
divergencies. For convenience, we regularize them by using conventional dimensional regu- 
larization (DREG). As is well known, DREG violates supersymmetry. In particular, SUSY- 
invariant physical quantities involving the gauge and Yukawa couplings g s and g s deviate by 
a finite amount in the exact SUSY limit. A SUSY-restoring counterterm can be derived from 
the Slavnov-Taylor identities This amounts to the following replacement of the Yukawa 
coupling g s in the gluino-quark-squark vertex 



37C. 



(4) 



The ultraviolet divergencies are removed by renormalization of the masses and the coupling. For 
the mass renormalization we use the on-shell scheme, while the strong coupling is renormalized 
in the MS scheme. After renormalization the virtual corrections still contain soft and collinear 
singularities which manifest themselves in double and single poles in £ = (d — 4)/2. They 
are cancelled by singular contributions from the four parton final states which are due to soft 
gluons and collinear massless partons. These contributions have to be computed analytically 
in d dimensions to perform the cancellation. For process ([Zb) we used the so-called dipole 
subtraction method [0], which has been recently generalized to include massive partons [@]. 
For process we used a variant of the phase space slicing method [£p. 

We now discuss our numerical results for the SUSY-QCD corrections, first for the process 
e + e~ — ► qqg. Fig. 4 shows the cross section c(E cut ) defined in Eq. (Q), where at NLO the 
unspecified part of the final state X can be a single gluon, two gluons or a massless qq pair. 
In all cases we require Ex > E cut and take in Fig. 4 E cut = 50 GeV. For the masses we choose 
nig = 400 GeV and mq = 300 GeV. The renormalization scale is set to /j = 1 TeV Fig. 4 shows 
that the SUSY-QCD corrections enhance the cross section in the peak region by about 20%. 

In Fig. 5a we show the dependence of o(E cut ) on the renormalization scale fj for fixed y/s = 
2 TeV and two different values for E cut . The /^-dependence is significantly reduced by the 
inclusion of the order oc^ corrections, thus leading to a more reliable theoretical prediction of 
the cross section. Fig. 5b shows the dependence of a(£ , cu t) on the cut parameter. We emphasize 
again that E cut is a physical parameter that can be chosen by the experimentalist. 

The inclusive cross section at LO and NLO for e + e~ — > qqg /qqg defined in (0) is shown in 
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Figure 4: Cross section for e + e — > ggX at leading order (dashed line) and including the SUSY- 
QCD corrections (full line) for m § = 400 GeV, m- q = 300 GeV, £ cut = 50 GeV and n = 1 TeV. 




Figure 5: In (a), the dependence of o(E cut ) on the renormalization scale is shown, where the 
reference scale is = 1 TeV. In (b), the dependence on E cut is depicted for /j = 1 TeV. The cm. 
energy was fixed to be y/s = 2 TeV, and masses were chosen as in Fig. 4. 



Fig. 6a and 6b for the cases niq < nig and niq > nig, respectively. The cross section o to t reaches 
values up to 5 fb for the case niq = 300 GeV and nig = 400 GeV. The SUSY-QCD corrections 
enhance the cross section in the peak region by about 25%. For niq > nig, the cross section is 
below 1 fb, and this case is probably of no phenomenological interest. The NLO corrections 
reach values up to 80%Q. The scale dependence of a to t is shown in Figs. 7a and 7b. For the 
phenomenologically more interesting case niq < nig the inclusion of the SUSY-QCD corrections 
significantly reduces the theoretical uncertainty due to the arbitrariness of the scale choice. 

'For niq > nig the process e + e~ — > qqg/qqg is only relevant for cm. energies between nig + nig and 2niq, since 
for yfs > 2m,q on-shell production of a squark-antisquark pair becomes possible. 
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Figure 7: Scale dependence of o to t in LO and NLO for the cases niq < nig (a) and niq > nig (b). 
The reference scale is /jq = niq + m g in both cases. 



4 SUSY-QCD corrections to e + e 
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Finally we study the impact of SUSY-QCD corrections on the production of 3 jets. For this we 
have to compute the SUSY-QCD loop corrections to the process e + e — > qqg. The 3-jet cross 
section is defined as follows: 



03 (y. 



cut; 



a(e + e ->3jets;yy>y cut J 



(5) 



where we take yy = (pt + Pj) 2 /s with p, denoting the four momenta of the final state partons. 
The cross section may be expanded as follows: 



a 3 (ycut) = ^c 3 + (-) 



.1,SUSY 1,QCD 



+ 0(a 3 s ). 



(6) 
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Figure 8: Impact of virtual SUSY-QCD corrections on the 3 -jet cross section. We require 
(Pi + Pj) 2 > (50 GeV) 2 for each pair (/, j) of the final state partons. 



Here C 3 ' D denotes the ordinary QCD correction at order a 2 . The 'genuine' SUSY-QCD 
corrections are defined by 



.l.SUSY 
2% o" 



Asvsy = ^^-q—. (7) 

The dependence of Asusy on Jcut turns out to be very small. We choose y cut = (50 GeM) 2 /s and 
plot in Fig. 8 Asusy as a function of y/s for niq = 300 GeV and nig = 400 GeV. Here we have 
used a decoupling scheme for the strong coupling, i.e. only the light quark flavours contribute to 
the running of a s . This ensures that Asusy — > as the squark and gluino masses become large 
compared to the cm. energy. Fig. 8 shows that the virtual SUSY-QCD corrections are tiny at 
the Z-pole, reach a maximum of about 1.5% around y/s = 1 TeV and eventually turn negative 
for large cm. energies. 



5 Conclusions 

The experimental verification of supersymmetric coupling relations is necessary to establish 
supersymmetry. The SU(3) coupling relation can be tested in e + e~ collisions by comparing 
measurements of the processes e + e~ — > qqg, qqg, qqg/qqg to precise, i.e. NLO calculations. 
We have computed the previously unknown order a 2 SUSY-QCD corrections to the above pro- 
cesses. The corrections enhance the LO cross sections for e + e~ — > qqg, qqg/qqg by about 20% 
or more in the peak region. The renormalization scale dependence is reduced for the cases of 
phenomenological interest. 
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